Systematic investigation of the elastic proton-deuteron differential cross section at 

intermediate energies 
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To investigate the importance of three-nucleon forces (3NF) systematically over a broad range of 
intermediate energies, the differential cross sections of elastic proton-deuteron scattering have been 
measured at proton bombarding energies of f08, f 20, f35, f50, f70 and f90 MeV at center-of-mass 
angles between 30° and f70°. Comparisons with Faddeev calculations show unambiguously the 
shortcomings of calculations employing only two-body forces and the necessity of including 3NF. 
They also show the limitations of the latest few-nucleon calculations at backward angles, especially 
at higher beam energies. Some of these discrepancies could be partially due to relativistic effects. 
Data at lowest energy are also compared with a recent calculation based on xPT. 
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During the last decade, the addition of three- 
nucleon forces (3NF) to modern high-quality nucleon- 
nucleon (NN) potentials, such as Nijmegen-I, Nijmegen- 
II, Reid93, CD-Bonn and Argonne-V18 (AV18) 1, 2, 3] 
has been shown necessary for many three-nucleon scat- 
tering observables, like the differential cross section and 
the vector analyzing power Aj, of elastic proton-deuteron 
scattering 0, H, H) S IS This necessity has been 
recognized before from the fact that 3N- and 4N- nu- 
clei are underbound by NN forces alone and also 
because low-lying spectra of li ght nuclei can not be de- 
scribed correctly without 3NF Il2l | . The most popu- 
lar 3NF are Urbana-i xjlij l ^ and a modified Tucson- 
Melbourne force TM' [Illllf (which no longer violates 
chiral symmetry). Various combinations of modern NN 
potentials and these 3NF have been worked out, fitted to 
the '^H binding energy ^^id subsequently apphed to 
3N scattering. The comparison to data revealed a mixed 
picture. In some cases NN forces alone work very well 
|l7| : in others, when NN force predictions fail, the inclu- 
sion of 3NF leads to a good or fairly good description 
0,0, 0,0, 0,0, and in stiU other cases neither NN forces 
alone nor the inclusion of these 3NF is sufficient to de- 
scribe the data H H H H lH . Recently, results of 
a systematic study of the nucleon vector-analyzing power 
of the reaction ^H(p,p(i) at several beam energies cover- 
ing a large range of center-of-mass angles were published 
[Tj . Serious discrepancies were observed at backward an- 
gles, showing that the spin structure of 3NF is not yet 
under control. Even though relatively precise data for 
the analyzing powers with large center-of-mass angular 
coverage exist in the literature 0, Q, 1 dif- 
ferential cross section of the reaction ^H()9, pd) few data 
sets are available. The existing data [iS 0, IU l23| . 



with the exception of the data from RIKEN |U, U| , are 
limited by low precision or small angular range. There- 
fore, in order to obtain a comprehensive picture of elas- 
tic proton-deuteron scattering at intermediate energies, a 
systematic measurement of the differential cross section 
of this reaction as a function of the center-of-mass scat- 
tering angle and the bombarding energy was performed 
at KVI. 

The present data are compared to state-of-the-art so- 
lutions of the Faddeev equations based on the above- 
mentioned high-precision NN potentials alone and in 
combination with the 3NF TM' and Urbana-IX. In ad- 
dition, we shall include very recent results using NN and 
3NF derived in chiral perturbation theory (xPT). This 
new approach to nuclear forces is an extension of on- 
going investigations in effective field theory applied to 
the nucleon itself, the ttN and the tt-tt systems. It is a 
systematic perturbative approach based on a smallness 
parameter, the ratio of generic external momenta to a 
certain mass scale of the order of the p mass. In the case 
of few-nucleon systems (and in the Weinberg scheme ^3) 
the nuclear forces are expanded in relation to that small- 
ness parameter (and pion-mass insertions). These forces 
are then inserted into the Schrodinger equation or equiva- 
lent formulations like those of Faddeev- Yakubovsky. Re- 
cently, nuclear forces up to next-to-next-to leading order 
(NNLO) have been worked out [11113. At this order, 
NN forces consist of one- and two-pion exchanges, which 
are parameter-free, and a string of contact forces of low 
chiral dimensions. The parameters of the latter (so-called 
low-energy constants) have been fixed by NN scattering 
data [13- With these parameters, the forces can describe 
the data up to about 200 MeV quite well. At this order 
3NF begin to arise |13| and consist of three different types 
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of topologies: a) a parameter-free two-pion exchange, b) 
a one-pion exchange between a nucleon and a contact 
force between the other two, and c) a pure three-nucleon 
contact force. As has been shown in |25j, there are only 
two parameters entering the 3NF of the types b) and c). 
They can be fixed by adjusting to the binding energy 
and to the doublet neutron-deuteron scattering length 
[2^. This new dynamical picture has already been suc- 
cessfully applied in 3N- and 4N- systems, especiall y a t 
lower energies for 3N scattering up to 65 MeV HHH^- 
It is now of great interest to see whether this approach 
will work even at 108 MeV, the lowest energy studied in 
this paper, or if higher orders in the chiral expansion or 
1/Mm corrections are needed. Other approaches to cal- 
culate effective three-body forces exist in the literature 
[iil E^ . At very low energies, techniques to solve the 
three-nucleon problem differently have been developed 
in recent years [s^, |^ • 

The experiment was performed at KVI using the Big- 
Bite Spectrometer BBS |23 in conjunction with the 
Euro-Supernova focal-plane detection system ESN [s^- 
Protons were obtained from either the KVI polarized ion- 
source js^l or the CUSP source and accelerated in the 
superconducting cyclotron AGOR js^- Measurements 
were made at bombarding energies of 108, 120, 135, 150, 
170 and 190 MeV. Since the lower acceleration limit of 
AGOR for polarized protons is at 120 MeV, protons with 
a kinetic energy of 108 MeV were obtained by passing 120 
MeV protons through an energy degrader. 

When polarized ions were used, the polarization was 
continuously measured using the KVI in-beam polarime- 
ter (IBP) in the high-energy beam-line. For these 
measurements, analyzing powers were obtained in addi- 
tion to cross sections. These analyzing powers were in 
total agreement with the published data j3| at the corre- 
sponding incident-beam energies. 

As targets, mixed solid CD2-CH2 targets with several 
thicknesses and a well-known CD2/CH2 ratio were used. 
The target thickness was determined from measurements 
of the differential cross section of elastic proton-proton 
scattering at several scattering angles. The measured 
differential cross sections of this reaction were compared 
with the results of NN calculations using Nijmegen-I, 
Nijmegen-II and Reid93 potentials, resulting in a normal- 
ization factor for the proton-deuteron differential cross 
sections. The analyzing power of the same reaction 
served as a cross-check of the polarization of the beam 
determined with IBP. 

The BBS was positioned at laboratory scattering an- 
gles between 17° and 50° in steps of 3°. At each scat- 
tering angle, deuterons and protons emerging from the 
reaction ^H(p,pd) were measured alternately. Deuterons 
were also measured at laboratory scattering angles of 5° , 
9°, 11° and 14°. In this way, center-of-mass angles be- 
tween 30° and 170° were covered for all six incident-beam 
energies. 



Measured differential cross sections are shown in figure 
n The statistical uncertainty, which is in general of the 
order of 2%, is smaller than the size of the symbols. The 
total systematic uncertainty, which is the quadratic sum 
of the uncertainty in the normalization factor, the uncer- 
tainty in the polarization where polarized protons were 
used, and the uncertainty in the detection efficiency, is 
in general ^ 7%, and should be considered as a point-to- 
point systematic uncertainty. 
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FIG. 1: Differential cross sections as a function of 6cm- 
The data measured in this work are plotted as open squares. 
The curves shown are calculations based solely on NN poten- 
tials (black band), calculations from Argonne-V18-|-Urbana- 
IX (solid line) and NN-I-TM' (gray band). At 135 MeV, data 
from 6] (circles), at 150 MeV data from i2gl (circles) and 
(diamonds) and at 190 MeV, data from (circles) and |23|| 
(solid squares) are shown as well. 



In all figures, the results of the Faddeev calculations 
using only two-nucleon interactions are shown as a black 
band. The width of this band represents the theoretical 
uncertainties in the calculations. Results from NN-I-TM' 
calculations are shown as a gray band. Further results 
from Argonne-V18-|-Urbana-IX are shown as solid lines. 
At 135 MeV, data from Sakai et al. Q are also shown. At 
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150 MeV, data from Postma et al. at 146 MeV |20| and 
Kuroda et al. at 155 MeV have been included in the 
figure. At 190 MeV, data from Adelberger and Brown 
[2^ at 198 MeV and Igo et al. at 181 MeV are 

shown. As can be seen, our data agree with most other 
existing data sets, taking into account the experimental 
uncertainties and the difference in incident-beam energy. 
At 135 MeV, our data lie systematically above the data 
from reference 0]. 

To make a better comparison between our data and 
the theoretical calculations shown in Fig. ^ the relative 
difference between them is shown in figure |21 In this 
figure, our data lie at zero and the calculations are shown 
by their relative deviations from our data. To avoid local 
fluctuations, these deviations are calculated by using a 
polynomial fit through the data points. 
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FIG. 2: Relative deviations of modern NN and NN-I-3N cal- 
culations from our data. The data points (open squares) lie 
at zero with the statistical uncertainties denoted at each data 
point. The total systematic uncertainty is shown as a light- 
gray band around zero. The meaning of the curves is the same 
as in figure Q 



In general, predictions based solely on NN interactions 
deviate from our data not more than w 60% over a range 



of three orders of magnitude. Inclusion of three-nucleon 
forces in the calculations leads to improved agreement at 
all incident-beam energies. With the high-precision data 
obtained in this work, which covers a broad center-of- 
mass angular region for six different incident-beam ener- 
gies, one can now study systematically the finer details 
of the calculations. 

At angles 9cm ~ 30°, almost all calculations overesti- 
mate the data. This is probably due to Coulomb effects 
not accounted for in the calculations. For 30° ^ 6'cm ^ 
60°, the predictions from different calculations with and 
without 3NF show some small disagreement among them- 
selves. Furthermore, the disagreement between the the- 
oretical predictions and our data is slightly outside the 
systematic uncertainties. At angles 9cm ^ 60°, the pre- 
dictions start to deviate from each other and from the 
data. In a large part of the angular range, especially at 
the higher bombarding energies, calculations using solely 
two-nucleon potentials fail to describe our data. This 
deviation is largest around 130° ^ ^cm ^ 150°. This 
angular range is part of the region of the minimum in 
the differential cross section and the place where three- 
nucleon force effects are expected to be largest 0j. The 
inclusion of three-nucleon forces in the calculations reme- 
dies these discrepancies, especially at the lower energies. 
However, starting at 135 MeV, deviations at backward 
angles around 6'cm ~ 140° can be observed. These de- 
viations increase with increasing the bombarding energy. 
Furthermore, a local minimum in the difference plot in 
figure [3 can be observed around 6'cm ~ 70°. This mini- 
mum is due to a shoulder that begins at 9cm ~ 60° and 
which is reproduced by the calculations but is more pro- 
nounced in our data. A similar pattern of disagreement 
exists in the vector- analyzing power 0. 

For the sake of clarity, we compare in Fig.|21our results 
for the lowest energy, 108 MeV, to predictions of chiral 
perturbation theory. The theory is shown as a band, 
which reflects the dependence on a momentum cut-off 
parameter underlying the specific effective field theory 
formulation used. The agreement with the data at this 
order of the theory is comparable to the results based 
on the conventional forces. Higher-order corrections are 
expected to improve the description of the data. 

The largest deviations of the theoretical predictions 
from our data occur at high incident-beam energies and 
large backward angles, i.e., a kinematic regime where the 
momentum transfer is largest. A plausible explanation 
for this is the higher-order effects, such as n-p or p-p 
exchange, which have not been included into the calcu- 
lations. Part of the disagreement may also be due to 
relativistic corrections which have not been properly 
taken into account in the calculations. 

In conclusion, the high-precision data presented in this 
work, covering a large center-of-mass angular region at 
several bombarding energies, make possible systematic 
study of the influence and the deflciencies of existing 
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FIG. 3: Absolute values and relative deviations of the cal- 
culated differential cross section in the framework of chiral 
perturbation theory (dark-gray band) in comparison with the 
data at 108 McV incident-beam energy. Also shown arc the 
data points from this work and, in the lower panel, the ex- 
perimental systematic uncertainty (light-gray band). 



three-nucleon forces. Calculations from NN-I-3N mod- 
els show deficiencies at backward angles at the higher 
bombarding energies. This may be an indication that 
the exchange of heavier mesons is missing in the calcula- 
tions. It may also be an indication that further relativis- 
tic corrections must be included. The calculations based 
on chiral perturbation theory for an incident-beam en- 
ergy of 108 MeV are promising and should be extended 
by incorporating N3L0 and 1/mjv corrections. 
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